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^ . ABSTRACT 

^ ' Photospheric stellar activity (i.e. dark spots or bright plages) might be an important source of noise and confusion in stellar radial- 

velocity (RV) measurements. Radial-velocimetry planet search surveys as well as follow-up of photometric transit surveys require a 
deeper understanding and characterization of the effects of stellar activities to differentiate them from planetary signals. We simulate 
dark spots on a rotating stellar photosphere. The variations in the photometry, RV, and spectral line shapes are characterized and 
analyzed according to the stellar inclination, the latitude, and the number of spots. We show that the anti-correlation between RV 
and bisector span, known to be a signature of activity, requires a good sampling to be resolved when there are several spots on the 
photosphere. The Lomb-Scargle periodograms of the RV variations induced by activity present power at the rotational period P,o, 
of the star and its two first harmonics -P„//2 and f ,-,„/3. Three adjusted sinusoids fixed at the fundamental period and its two-first 
(— I ' harmonics allow us to remove about 90% of the RV jitter amplitude. We apply and validate our approach on four known active planet- 

^ ^, host stars: HD 189733, GJ674, CoRoT-7, and iHor. We succeed in fitting simultaneously activity and planetary signals on GJ674 

I . and CoRoT-7. This simultaneous modeling of the activity and planetary parameters leads to slightly higher masses of CoRoT-7b and 

O ' c of respectively, 5.7 ± 2.5 M£„„/i and 13.1 ± 4.1 M£„„/i- The larger uncertainties properly take into account the stellar active jitter. 

We exclude short-period low-mass exoplanets around t Hor. For data with realistic time-sampling and white Gaussian noise, we use 
C/3 , simulations to show that our approach is effective in distinguishing reflex-motion due to a planetary companion and stellar-activity- 

■ induced RV variations provided that 1) the planetary orbital period is not close to that of the stellar rotation or one of its two first 

harmonics, 2) the semi-amplitude of the planet exceeds ~30% of the semi-amplitude of the active signal, 3) the rotational period of 
the star is accurately known, and 4) the data cover more than one stellar rotational period. 

^ ' Key words, techniques: radial velocities - stars: activity - stars: individual: t Hor, HD 189733, GJ674, CoRoT-7 
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1. Introduction The anti-correlation that exists between RV and bisector 

span (the measurement of the spectral line asymmetry) is 

, . . ,. , , . , ., , , currently used as a signature of RV variations induced by 

Hish-precision radial-velocimetry has untd now been the most . „ /~v i » i inm u » » i oaaon 
■ ^ ^ ,. . , stellar- activity (e.g. Oueloz et al. 2001; Hatzes et al. 2002). 

erhcient way or discovering planetary systems. However, an „ • .u- ? w » u» » »u • j j 

. ,•: , • , , , , , , , • , Removing this anti-correlation to subtract the activity-induced 

active star displays on Its photosphere dark spots and/or bright „ » j • » / e 

^ ^ . , , fr,, . , ^ . . ° RV may allow us to derive more accurate parameters or the 

plages that rotate with the star. These inhomogeneities across , . . . , . , . , onm d • * i oaahn 

^ ° „ ^ . , ° , . , planetary system (e.g. Melo et al. 2007; Boisse et al. 2009). 

the stellar surface can induce RV shirts due to changes in the a »u u ,. ^- ,.■ ■ u- * n j i 

^ . ■ , , , Another approach to distinguishing Stellar activity and planetary 

spectral line shape that may add contusions with the Doppler • . ■ ^ e j- \.- » t/^ i • • i X 

, ^„ . , , . . ^ , , signatures consists or adjusting two Keplerian signals, one with 

H reflex-motion due to a planetary companion (e.g. Queloz et al. , . ■ j j »u j \.u »u » n » i 

d ^o/M • . . ^r^T^A . ^nr-c-Trr,, ,• , the planetary period and the second with the stellar rotational 

2001; Desidera et al. 2004; Huelamo et al. 2008). The amplitude ■ ■/ n ,.1 oaat c n » 1 oaaon 

^, ' , , , , , . . ^, , ^ , period (e.g. Bonfils et al. 2007; Forveille et al. 2008). 
or the RV shirts depend on the v sin ; or the star, the spectrograph 

resolution and both the size and temperature of spot (Saar & 

Donahue, 1997; Hatzes 1999; Desort et al. 2007). Stellar activity 

generates RV jitter and prevents accurate measurements of the In these past few years, the impact of the stellar activity 

stellar motions. For these reasons, active stars are then usually has been studied. Desort et al. (2007) quantified how the 

discarded from RV surveys using criteria based on the activity impact of stellar spots on RV measurements depended on the 

index R^j^ and/or vsin/. However, photometric transit-search spectral type, the rotational velocity, and the characteristics 

missions (such as CoRoT and Kepler) require RV measurements of the spectrograph used. Santos et al. (2009) investigated 

to establish the planetary nature of the transiting candidates long-term HARPS measurements of a sample of stars with 

and characterize their true masses. These candidates include known activity cycles. Lagrange et al. (2010) and Meunier et 

active stars adding strong confusions and difficulties in the RV al. (2010) analyzed spots, plages, and convection simulations 

follow-up. computed using solar data. 
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In this paper, we characterize variations in all the observed 
parameters derived from RV measurements as well as from pho- 
tometric measurements due to stellar activity. We study how the 
Keplerian fit used to search for planets in RV data is affected 
by spots and we test an approach to subtract RV jitter based 
on harmonic decompositions of the star rotation. For this, we 
use simulations of spectroscopic measurements of rotating spot- 
ted stars and validate our approach on active stars monitored by 
high-precision spectrographs: HD 189733 (Boisse et al. 2009), 
GJ 674 (Bonfils et al. 2007), CoRoT-7 (Queloz et al. 2009), and 
iHor(Vauclairet al. 2008). 

2. Simulations of activity-induced radial velocity 

2.1. SOAP tool: dark spot simulations 

SOAP (Boisse et al., in prep.) is a program that calculates the 
photometric, RV, and line shape modulations induced by one (or 
more) cool spots on a rotating stellar surface. SOAP computes 
the rotational broadening of a spectral line by sampling the stel- 
lar disk on a grid. For each grid cell, a Gaussian function rep- 
resents the typical line of the emergent spectrum. The Gaussian 
is Doppler-shifted according to the projected rotational velocity 
(v sin /) and weighted by a linear limb-darkening law. The stel- 
lar spectrum output by the program is the sum of all contribu- 
tions from all grid cells. The spot is considered as a dark surface 
without any emission of light, so we cannot compute different 
temperatures for the spot. For a given spot (defined by its lati- 
tude, longitude, and size), SOAP computes which grid cells are 
obscured and removes their contribution to the integrated stellar 
spectrum. 

For the simulation, we choose a GOV star with a radius of 
1.1 Rq and a vsin/=5.7 kms"', a linear coefficient of the limb 
darkening of 0.6, and a spectrograph resolution of 110 000 in 
order to be in the same conditions of the l Hor data presented in 
Sect. 3.4. We fixed arbitrarily a dark spot size of 1% of the vis- 
ible stellar surface. The stellar spectrum output by SOAP is an 
averaged spectral line, equivalent to the cross-correlation func- 
tion (CCF) computed to measure RV with real data. We fit the 
simulated CCF with a Gaussian that sets the parameters of the 
CCF (RV, contrast, FWHM, and bisector span). The photometric 
flux is also computed. 

2.2. RV variations due to a dark spot 

Fig. [U shows the RV modulations due to a spot as a function of 
time for different inclinations / of the star with the line of sight 
and different spot latitudes lat. These two parameters clearly 



modify the pattern of the RV modulation. If the spot remains vis- 
ible during all the stellar rotation (lat > i), the shape is close to 
a sinusoidal function (Fig.[Tl right). If the spot is hidden during 
the rotation of the star (Fig.[Tl left), the RV variation resembles 
a Rossiter-McLaughlin (RM) effect (Rossiter 1924; McLaughlin 
1924). 

Fig. |2] shows the Lomb-Scargle periodograms of the three 
cases showed in Fig. [T] Main peaks are clearly detected at the 
rotational period of the star Pror, as well as the two first har- 
monics Protl'2' and Protl^- We note that the energy in each peak 
varies with the shape of the RV modulation. Multiples of the 
rotational period are never found. Low-amplitude signal is de- 
tected at Proi/4- but only when the star is seen equator-on and the 
spot is close to the equator. In that case, the RV change departs 
strongly from a sinusoidal shape and the periodogram exhibits 
a stronger amplitude excess at Proi/2 rather than at the stellar 
rotational period. 
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Fig. 2. Lomb-Scargle periodograms of the three RV modulations 
showed in Fig. [1] The fundamental frequency, Proi, and its first 
harmonics are detected. 
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Fig. 1. RV modulations due to one spot as a function of time 
(expressed in rotational period unit). At f =0, the dark spot of 1 % 
of the visible stellar surface is in front of the line of sight. The 
shape of the signal changes with the inclination / of the star and 
the latitude lat of the spot, labelled in the top left of each panel. 



A third of the active regions, where spots grow and decay, 
appear at the same location as a previous active region and their 
lifetime can be several rotation timescales (e.g. Howard 1996). 
Hence, the phase of the RV jitter is preserved when the spot 
movement is only linked to the stellar rotation. In Fig. |3] we 
simulated the RV modulation taking into account the evolution 
of a spot, i.e. when the spot size and/or temperature changes with 
time. The Lomb-Scargle periodogram has identical peaks at Prot 
and its two first harmonics. Finally, the periods detected in the 
periodogram are the same for the following configurations: 1) a 
star with different inclinations, 2) spots at different latitudes, 3) 
spot size varying with time, and 4) several spots on the stellar 
surface (cf. Sect. 
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Fig. 3. Top: RV as a function of time when the spot size de- 
creases. Bottom: Lomb-Scargle periodogram of the RV. Peaks 
correspond to the rotational period and its harmonics. 



2.3. CCF parameter variation due to a dark spot 

The bisector span (BIS) is a measurement of the asymmetry of 
the CCF, which corresponds more or less to the average line of 
the spectrum. An anti-correlation between the RV and the BIS is 
a signature of activity-induced RV variations. The slope of the 
anti-correlation depends on the spot size, the v sin / of the star 
and the resolution of the instrument (Santos et al. 2003; Desort 
et al. 2007). 

The BIS was defined by Queloz et al. (2001) by computing 
the RV difference of the upper and the lower part of the bisector. 
This measurement appears unreliable at low signal-to-noise ratio 
(SNR). We describe here another approach to estimate the asym- 
metry of a CCF. We fit the CCF in two steps selecting the upper 
and bottom part to compute the RV;„g/, and RV;„„., respectively. 
Where cr is the width of the CCF, the upper part of the CCF is 
defined in the range [-oo:-lcr][-i-lcr:-i-oo] and the lower part is de- 
fined in the range given by [-oo:-3cr][-lcr:-i-lcr][H-3cr:H-oo]. The 
RV/ow is more sensitive to the variation in the bottom of the line, 
and then, more sensitive to stellar activity as shown in Fig. |4] 
The difference RV -RV/oh, - ^ span gives a measurement of 
the CCF asymmetry. In Fig.|5] the two techniques of CCF asym- 
metry measurement are plotted for two regimes of SNR illustrat- 
ing that at low SNR, the V span is more robust than BIS. 

For one spot, the V span or BIS is anti-correlated with the 
RV as expected (see Fig. |6]l. Although the shape of the anti- 
correlation is not a straight line, but a loop more or less like an 
inclined eight. Depending on the geometrical configuration, the 
mean slope and the width of the eight-\\ke loop vary. 

The subtraction of the average slope of the anti-correlation 
has previously been used to remove the activity-induced RV in 
order to improve the parameters of a Keplerian planetary orbit 
(e.g. Melo et al. 2007; Boisse et al. 2009). We simulated the ef- 
fect of this approach and found that the pattern of the RV resid- 
uals after the correction is identical with a lower amplitude of 
the initial RVs. The detected periods in the Lomb-Scargle peri- 
odogram of the RV residuals are indeed the same, with an am- 
plitude ratio equivalent to those in the initial RVs periodogram. 
The use of the RV-BIS anti-correlation to correct RV jitter has 
several limitations. First, it merely reduces the amplitude of the 
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Fig. 4. RV/jjj/, (black) and RV/o,v (red) as a function of time. The 
star has an inclination /=40° and the latitude of the spot is equal 
to +30°. RV/oit. is more sensitive to activity than RV^/^;,. 



Vspan/BIS [m/s] 
40 



BIS [m/s] 





1 1 1 1 1 1 1 _ 




-0.41 










- 1 , 


A - 

, 1 , , , , 1 - 


-50 


50 



100 

50 

-50 
-100 



L^.'-A. .... 


1 1 1 . 




0.12 


■ — 






, , 1 : 


-50 


50 




RV [m/s] 



RV [m/s] 



-50 
RV [m/s] 



Fig. 5. Left: BIS (triangles) and V span (squares) as a function of 
RV for a simulated spot rotating at latitude -1-50° on a star with an 
inclination /=40°. Middle and Right: Respectively, BIS (Queloz 
et al. 2001) and V span (this paper) as a function of RV for the 
same simulated spot with 20ms^' additional photonic noise in 
the CCF. The lines are the least squares fit to the data. The num- 
bers in the right-hand corner is the value of the slope of the fit. 



jitter signal, and does not alter its frequency content. Second, it 
can only be used if the vsini of the star is sufficient to resolve 
the effect of spots on the bisector (Desort et al. 2007; Bonfils et 
al. 2007). Third, we have shown that the degree of coiTelation 
depends on the spot latitude and stellar inclination. Finally, the 
coiTelation may be absent if effects other than activity contribute 
significantly to the RV variations (Queloz et al. 2009). 

Our simulation does not take into account the variation in 
line width with stellar effective temperature, which is affected 
by the presence of spots as measured and calibrated by Santos 
et al. (2002). In our simulated CCF, the flux contribution of the 
area occulted by the spot is subtracted from the mean stellar line. 
The resulted CCF is distorted relative to a Gaussian profile. As 
the spot rotates with the stellar surface, the contribution of the 
spot moves inside the CCF, from the blue to the red wings. The 
contrast and FWHM of the fitted CCF vary in anti-coiTelation as 
illustrated in Fig. |7]where the CCF parameters are expressed as 
a percentage of variation relative to the mean value. When the 
spot is in front of the line of sight, its contribution affects the 
center of the line profile, the contrast is reduced and the FWHM 
is larger. In contrast, when the spot contribution is in the wings 
of the line, the FWHM is smaller. In practice, CCF contrast may 
be more sensitive to both the conditions of observation (e.g. sky 
background) and the diffuse light in the spectrograph. The CCF 
FWHM is then probably the most reliable parameter to charac- 
terize an effect due to the activity. 
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Fig. 6. ^ span as a function of RV for the three RV modulations 
of Fig.l. The relation is not a perfect anti-correlation but looks 
more like an inclined eight shape. 
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Fig. 7. From top to bottom: BIS, FWHM, contrast of the CFF 
and the photometric flux as a function of time. The variation 
in the parameters of the CCF are expressed as a percentage of 
the variation around the mean value. The star has an inclination 
i-4Q° and the latitude of the spot is equal to +30° (RV shown in 
Fig.\l\middle). Att-l, the spot is in front of the line of sight. 

2.4. Correlation witli piiotometry variations 

We compared the variations in the photometric flux with those of 
the CCF parameters (asymmetry, width, and contrast) for differ- 
ent stellar inclinations and spot latitudes (see Fig. [T). The con- 
trast is anti-correlated with the FWHM. Two regimes are ob- 
served between the CCF contrast and the photometric flux (see 
Fig.[8]l. The simulations compute the distortion of the CCF as a 
lack of flux masked by the spot. When the spot is in front of the 
line of sight, i.e. when the emitted flux is at its lowest, there is 
a bump at the line center and the wings appear broader. In con- 
trast, when the spot is in the limb, i.e. when the distortion is in the 
wings and the emitted flux is at its medium value, the line width 
appears to be naiTower. This explains why an anti-coiTelation is 
observed between CCF width and photometry. The other part, 
where a correlation is observed between the parameters is not 
well understood. As shown in Fig. |7] the contrast and FWHM 
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Fig. 8. FWHM of the CCF as a function of the photometric flux 
for the three cases of Fig. [1] The variations in the FWHM are 
expressed as a percentage of the variation around the mean value 
(here, equal to zero). 

vary in anti-correlation, the effect in terms of FWHM then also 
depends on the contrast variations. This part of the relation be- 
tween FWHM and photometry might originate from our inabil- 
ity to compute the different temperatures of the spot. On Corot-7 
(Queloz et al. 2009), the photometric measurements from the 
Swiss Euler telescope are anti-coiTelated with the HARPS CCF 
width. 

2.5. RV fit of a dark spot 

We now attempt to remove or at least reduce the stellar activity 
signals in order to identify a planetary signal hidden in the RV 
jitter. We add to the simulated RV time series a Gaussian noise of 
0--1 ms"' to take into account photon noise and/or instrumental 
noise. The main shape of the RV modulations and the peaks in 
the periodograms are not affected with such a noise having an 
amplitude 100 times lower than the activity jitter. 

The Lomb-Scargle periodogram corresponds to sinusoidal 
decompositions of the data. We saw previously that the main 
signal are at the stellar rotational period Proi and its two first 
harmonics Protl'2' and Protl^- In some configurations, signal is 
present at Protl^- We firstly fit five sinusoids to the three simu- 
lated cases of Fig. [1] with periods fixed at the rotational period 
Prnt, and its first harmonics. The residuals are plotted in Fig. |9] 
and the Lomb-Scargle periodogram of the residuals in Fig. [10] 
The parameters of the sinusoids fitted are reported in Table 1. 
The active jitter is reduced from more than 100 ms"' peak-to- 
peak to less than 3 ms"'. In reality, only three sinusoids are 
needed to decrease the semi-amplitude of the RV jitter by more 
than 87% because the semi-amplitude of the fourth (P,o//4) and 
fifth {Proi/5) sinusoids are negligible in almost all cases. 

We compared the fit of activity signals with Keplerian and 
sinusoids. Fitting a Keplerian is equivalent to fitting a sinusoid at 
the same period and its first harmonics. The parameters of the fit 
are reported in Table 2. The number of free parameters is equal to 
16 for a fit with three Keplerians and to 10 with three sinusoids. 
Taking into account only the fundamental period Pr„i as a free 
parameter (and not the other periods fixed as harmonics), there 
are 14 free parameters for a fit with three Keplerians against 8 
for a fit with three sinusoids. Moreover, if some energy due to 
activity remains in the residuals, it is at the next harmonic for a 
fit with sinusoids, while it may be at any other harmonic for a 
fit with Keplerians. This is why we choose to fit the active jitter 
with three sinusoids fixed at the stellar rotational period and its 
two first harmonics. 

The latest harmonics may appear because the sampling of 
the simulation is almost perfect. In the following section, we 
checked that in practice when the sampling and coverage of 
the stellar activity signal are limited, a three sinusoid fit with 
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fixed periods at the rotation period, and its two first harmonics 
is sufficient. 
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Fig. 9. Residuals from the five-sinusoid fit as a function of time 
of the RV variations showed in Fig.[T] 




Fig. 10. Lomb-Scargle periodogram of the residuals of the 3 RV 
variations showed in Fig. [T] after the five-sinusoid fit with fixed 
periods at the rotational period of the star and its harmonics. 



2.6. Several dark spots 

We simulated and studied the apparent shifts in RV induced by 
several spots on the stellar surface. Two spots of 1 % of the stellar 
surface are successively shifted by 30, 60, and 120° in longitude 
on an inclined star with / = 40°. One spot is fixed at latitude H-30° 
and the other at the equator. Differential rotation is not consid- 
ered. The three RV signatures and periodograms are plotted in 
Figs. [TT] and [121 Only the three periods Pro,/2, and Protl'i 
are identified. Three sinusoids fixed at these periods fit the RV 
modulation and reduce the amplitude of the jitter by up to 90%. 

For our three simulated cases with two spots, the anti- 
correlation between Nspan and RV is more or less observed on 
Fig. [13] but may be affected by one spot possibly compensat- 
ing for the other. We emphasized that the slope of the anti- 
correlation is shallow and requires data with a high sampling 
to be resolved (at least 10 measurements per rotation period). 
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Fig. 11. RV as a function of time induced by two spots on the 
stellar surface. One spot is at latitude -1-30°, the other at the equa- 
tor. From top to bottom: the shape of the variation varies with the 
difference in longitude between the two spots, labelled in the top 
left of each panel. 



The pattern observed for one spot between the CCF FWHM 
and the photometry can be strongly modified by the presence 
of several spots as shown in Fig. [14] It may not display a clear 
signature and the anti-correlation may be barely seen. However, 
as both parameters are not modified by the presence of a planet 
(except for a transit), this anti-correlation due to stellar activity 
would always be observable. This is in contrast to the RV and 
^ span (or BIS) anti-correlation, which may be hidden by a planet 
in the RV data. 



2.7. Bright spots 

We consider only dark spots, but other inhomogeneities are 
present on the surface of an active star. For instance, plages are 
bright structures that may be related to dark spots in active re- 
gions. We simulated bright spots and characterized the induced 
variations in the RV, CCF parameters, and photometry. We found 
that the amplitudes of the variations are identical for character- 
istics similar to those of a dark spot (except the brightness), but 
the shapes of all parameters are reversed. The Vspan and the RV 
then remain anti-correlated. Similarly, the relation between pa- 
rameters of the CCF and the photometry persists. 



3. Application to real data 

3.1. HD 189733 

The active K2V star HD 189733 and its transiting planetary com- 
panion, which has a 2.2-day orbital period, was monitored by 
Boisse et al. (2009). They used the high-resolution spectrograph 
SOPHIE mounted on the 1 .93 -m telescope at the Observatoire de 
Haute-Provence to obtain 55 spectra of HD 189733 over nearly 
two months. The RV measurements subtracted from their fit of 
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Fig. 12. Lomb-Scargle periodograms of the three RV variations 
showed in Fig. [TT] when two spots are on the stellar surface. Only 
the rotational period of the star and its two-first harmonics are 
detected. As noted in Sect. 12.21 the maximal peak is not at Pr„i /2, 
so there is not a significant detection of the harmonic Pnnl^- 
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Fig. 13. span as a function of RV for the three RV variations 
showed in Fig.[TT]when two spots are on the stellar surface. The 
square points represent a sampling with one measurement each 
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the planetary companion are variable because of the activity of 
the star. HD 189733 has a stellar rotational period of 11.953d 
(Henry & Winn 2009). We computed the Lomb-Scargle peri- 
odogram of the residuals from the Keplerian fit in the top panel 
of Fig. [15] The fundamental period at P,.o, and its two first har- 
monics (Pra,/2 and Praf/3) are detected with false alarm prob- 
ability lower than 10"'. False alarm probabilities are computed 
by performing random permutations of the data, keeping the ob- 
serving time fixed, as described in Lovis et al. (2010). There is 
no significant signal at other harmonics. The dispersion value 
of 9.1 ms ' are reduced to 5.5 ms ' when the residuals are fit- 
ted by three sinusoids fixed at the given periods. This dispersion 
is comparable to that obtained by the RV-BIS anti-correlation 
correction applied by Boisse et al. (2009). The main limitation 
here is the intrinsic precision of the spectrograph. The Lomb- 
Scargle periodogram of the residuals after the simultaneous fit of 
a Keplerian and three sinusoids with periods fixed at Prot and its 
two first harmonics is shown as a dashed line at the top of Fig.fTSl 
It shows no evidence of another companion in the system taking 
into account the current SOPHIE accuracy (^ 4-5ms '). 

We show the effect of the sampling in the middle and bottom 
panels of the Fig. [15] We simulated a dark spot on the star with 
the closest parameters to those of HD 189733. Spot parameters 
(0.45% of the visible stellar surface and lat—+liO°) were selected 
to ensure that the simulated photometric variations are compa- 
rable to those observed by MOST (Boisse et al. 2009). In the 
bottom panel, the Lomb-Scargle is computed for a quasi-perfect 
sampling during which a measurement is taken each Pro//1000. 
In the middle panel, the sampling is that of the measurements 
of HD 189733 in Boisse et al. (2009). The pattern of the Lomb- 
Scargle periodogram for the simulations with the sampling of 
the SOPHIE data {middle panel) is comparable to those of real 
data (top panel). As a dashed line, we superimpose the Lomb- 
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Fig. 14. FWHM of the fitted Gaussian as a function of the pho- 
tometric flux of the three RV variations shown in Fig. [TT] when 
two spots are on the stellar surface. The square points represent 
a sampling with one measurement each Pror/lO. 
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Fig. 15. Top: Lomb-Scargle periodogram of the residuals from 
the Keplerian fit of the SOPHIE RV of HD 189733 (black curve). 
The three peaks at the rotational period and the two first harmon- 
ics are detected. The Lomb-Scargle periodogram of the residu- 
als after the simultaneous Keplerian fit and the harmonic filter- 
ing of the stellar activity is indicated by a dashed line. Middle: 
Lomb-Scargle periodogram of a simulation of one dark spot 
on the surface of a star with parameters closest to those of 
HD 189733 (black curve). The sampling is identical to that of the 
SOPHIE data for HD 189733. The lomb-Scargle periodogram of 
the residuals after the simultaneous Keplerian fit and the har- 
monic filtering of the stellar activity is indicated by a dashed 
line. Bottom: Lomb-Scargle periodogram of the same simula- 
tion as in the middle panel with a very good data sampling (one 
measurement each Pror/1000). 



weaker dispersion in the residuals, cr((9 - C)=0.65 ms"' instead 
of cr(0 - C)=0.82 ms which is closest to the cuiTent HARPS 
accuracy and equal to the uncertainty in each measurement, and 
we were able to reduce the;^'^ to 1.36. 

For the original RV data, we computed the Lomb-Scargle 
periodogram, (Fig. [17] top), which shows that the main period 
is that of the planet. After removing the RV variations due to 
the planet, the periodogram of the residuals displays periodici- 
ties at the rotational period and its two first harmonics (Fig. [TT] 
bottom). After the simultaneous fit of a Keplerian and the har- 
monic filtering, the periodogram shows neither a residual nor an 
enhanced periodicity. 
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Fig. 16. Fit of the HARPS RV of GJ674 with one Keplerian 
and three sinusoids fixed to the rotational period (top 
right P,o,=34.85d) and the two first harmonics (bottom 
left P„„/2=17.425d and bottom right P„,/3=l 1.6167d). The 
Keplerian fit (top left) gives the planetary parameters. Each plot 
is given as a function of orbital phase after removing the other 
three signals. Individual error bars are also plotted. 



Scargle periodogram of the residuals of the harmonic filtering 
of the simulated RV with the SOPHIE data sampling. No pe- 
riodicities are found to have greater amplitudes because of the 
harmonic filtering beyond the noise level. 

3.2. GJ674 

GJ 674 is a moderately active M2.5 V dwarf hosting a planet 
with a 4.69-day period (Bonfils et al. 2007). A superimposed 
signal with a periodicity of roughly 35 days is also visible in 
the RV measurements. Bonfils et al. (2007) demonstrated that 
their variations coincide with the stellar rotation period by ana- 
lyzing photometry and active lines (Ca II and H a). These signals 
originate in active regions rotating with the stellar surface. They 
described the data with a two-Keplerian model with a 4.69-day 
period planet and a Keplerian to fit the active feature. In Fig.[T6l 
we fit the RV derived from HARPS spectra with three sinusoids 
each of period fixed to be the rotational period (Pro,=34.85d) and 
the two first harmonics (P„,/2=17.425d and P„«/3=l 1.6167d) 
and one Keplerian that provides the planetary parameters. These 
fitting results are in agreement with Bonfils et al. (2007) (cf. 
Table 3), being no significant change in either the fitted pa- 
rameters or the error bars. On the other hand, we obtained a 



Table 2. Keplerian orbital solutions for GJ674b planet as pub- 
lished in Bonfils et al. (2007) and fitted in this paper. 



Planetary parameters Bonfils et al. (2007) This paper 



Pp [days] 


4.694 + 0.007 


4.694 ± 0.002 


K[ms-'] 


8.70 ±0.19 


8.9 ± 0.3 


e 


0.20 ± 0.02 


0.19 ±0.03 


oj [deg] 


143 ±6 


159 ± 10 


To [JD] 


53780.09 ± 0.08 


53780.25 ± 0.12 


)!J2sini'' [Me] 


11.09 


11.39 


o-^o-c)'' [ms"'] 


0.82 


0.65 


reduced x~ 


2.57 


1.36 



" assuming = 0.35 Mq (Bonfils et al. 2007) 
(T^o-c) after the fit 



3.3. CoRoT-7 

The photometric transit search with the CoRoT satellite has re- 
ported the discovery of a planetary companion CoRoT-7b around 
an active V= 11.7 G9 V star (Leger et al. 2009) with an orbital pe- 
riod of 0.85 days. Queloz et al. (2009) reported on the intensive 
campaign carried out with HARPS at the 3.6-m telescope at La 
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Fig. 17. Top: Lomb-Scargle periodogram of the RV of GJ674. 
The main peaks are due to the periodicity of the planet and its 
one-day alias. Bottom: Lomb-Scargle periodogram of the resid- 
uals from the Keplerian fit of the GJ674 RV data (black curve). 
The main peak is the rotational period of the star. The two first 
harmonics are also visible. The Lomb-Scargle periodogram of 
the residuals after the simultaneous Keplerian fit and harmonic 
filtering of the stellar activity is shown as a dashed line. 



which were fixed to be zero. The differences from the published 
values of the periods are smaller 0.5% and of the transit phase 
of CoRoT-7b is less than 0.2% of the transit period. To measure 
the semi-amplitude and then the mass of the planets, we fixed 
the period and the TO of the transiting companion. We then fit- 
ted simultaneously three sinusoids for the active jitter and two 
Keplerians with null-eccentricity. We obtained the parameters 
given in Table[3] The period of 3.70 ± 0.02 d found for CoRoT- 
7c agrees with the value of 3.698 + 0.003 d found by Queloz et 
al. (2009). The residuals of 3.1 ms ' and the reduced which 
is equal to 2., are comparable to those obtained by Queloz et 
al. (2009) considering that we have used fewer data. We note, 
however, that the semi-amplitude of the planets 4.5 ± 0.7 ms"' 
for CoRoT-7b and 6.1 + 0.6 ms ' for CoRoT-7c are larger than 
in Queloz et al. (2009). For comparison, the same study is per- 
formed on another data set of bjd=[54775. 8-54807. 3], with a 
total of 34 measurements. The two planets are detected in the 
Lomb-Scargle periodogram with false alarm probabilities lower 
than 10"^, whereas no other periods are detected with false alarm 
probability lower than 0.5 (except for the one-day alias due to 
the data sampling). The differences from the published values 
are 4.5% for the CoRoT-7c period, 2% for the CoRoT-7b tran- 
sit period, and 9% for the CoRoT-7b transit phase. The semi- 
amplitude for the planets are 2.9 +1.1 ms ' for CoRoT-7b and 
4.6 + 0.9 ms-' for CoRoT-7c. 



Silla. The RV variations are found to be dominated by the activ- 
ity of the star, i.e. rotational modulation from cool spots on the 
stellar surface with an estimated period close to 23 days. Two ap- 
proaches were used to disentangle the Doppler motion from the 
active jitter. They each detected two signals: the CoRoT-7b tran- 
sit period and a second planet CoRoT-7c at a period of 3.69 days. 
The authors inferred a mass of 4.8 + 0.8 M^a,,/, for CoRoT-7b 
and, assuming that both planets are on coplanar orbits, a mass 
of 8.4 + 0.9 M.Earth for CoRoT-7c. The second of the two ap- 
proaches used in Queloz et al. (2009) is a modeling of the active 
jitter by a harmonic decomposition of the rotational period. The 
authors subtracted their model from the RV data before detecting 
and characterizing the planetary system. Here, we wish to fit si- 
multaneously the effect of both activity and the planetary system 
on the data as we did previously for GJ674. 

We used the 37 last days of HARPS data when the sampling 
had a higher frequency on a long timescale (several points each 
night and a total of 59 measurements) and in order that the dis- 
tribution of spots on the stellar surface did not change too much 
(and therefore the parameters of the fit for the active signal). 
This assumption is consistent with that of the photospheric ac- 
tivity study of Lanza et al. (2010) for the CoRoT photometric 
measurements. They found that active regions occur on three 
stable active longitudes, whose overall lifetime may exceed the 
duration of the light curve. They note that the large active re- 
gions evolve on timescales ranging from two weeks to a few 
months. However, we show in Sect. 2 that an evolving active 
region is also closely fitted by our harmonic filtering. We fitted 
simultaneously three sinusoids to the active jitter with periods 
fixed at the rotational periods (23d.) and its two first harmonics 
(11. 5d and 7.66667d.). The Lomb-Scargle of the residuals shows 
a clear peak near 3.69 d and another one near 0.85 d with false 
alarm probabilities lower than 5.10""* (Fig.fTSl). No other signif- 
icant periods are detected with false alarm probability greater 
than 0.5. We fitted simultaneously three sinusoids to the active 
jitter and two Keplerians to the possible companions. No pa- 
rameters were fixed for the Keplerians except the eccentricities. 
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Fig. 18. Lomb-Scargle periodograms of the 37 last days of 
HARPS RV data of CoRoT-7. The dotted line shows the level of 
the mean uncertainty in RV measurements, a) The main peak in 
the raw data is owing to the rotational period of the star (~ 23c/.). 
b) After the harmonic filtering of the data, the highest peak of 
the residuals corresponds to the CoRoT-7c period, c) The resid- 
uals of the simultaneous fit of the activity and one Keplerian for 
CoRoT-7c show two peaks, the highest at the CoRoT-7b period 
and the other at its one-day alias, d) No significant periods re- 
mains above the level of the noise after the simultaneous fit of 
the activity and of two Keplerians for the possible companions. 
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Table 4. Adopted values for the two CoRoT-7 planets as derived 
from this paper. 



Parameters 


CoRoT-7b 


CoRoT-7c 


Pp [days] 


0.8536 (fixed) 


3.697 ± 0.019 


K[ms-i] 


4.0 ± 1.6+ 


5.6 ± 1.6"'" 


e 


(fixed) 


(fixed) 


To [JD] 


54899.761 (fixed) 


54899.2 ± 0.7 


mi [Me] 


5.7 + 2.5* 


13.2 + 4.1* 



+ Error bars in semi-amplitudes take into account a systematic noise 
due to active jitter of 1.5 ms"'. 

* Error bars in masses take into account the stellar mass uncertainty 
and a systematic noise due to active jitter. 



Our method is robust, but these differences illustrate the dif- 
ficulty in measuring the amplitudes accurately in the presence of 
activity. Moreover, it led us to conclude that the error bars for 
the semi-amplitudes are underestimated in our study as in that 
of Queloz et al. (2009). We infer approximately that a system- 
atic noise due to active jitter of 1.5 ms"' must be added quadrat- 
ically to the error bars. We then determine a semi-amplitude of 
4.0+ 1.6ms-i for CoRoT-7b and 5.6+ 1.6 ms"' for CoRoT-7c. 
The masses are 5.7 ±2.5 MEarth for CoRoT-7b, which agrees 
with the value of Queloz et al. (2009), and 13.2 + 4.1 MEarth for 
CoRoT-7c (assuming that both planets are on coplanar orbits), 
slightly higher than its published value. The error bars then ac- 
count for a systematic noise due to active jitter and the stellar 
mass uncertainty (0.93 + 0.03 Mq). By evaluating a weighted av- 
erage, we find for CoRoT-7c a period of 3.697 + 0.019 day, in 
agreement with Queloz et al. (2009). These new values for the 
CoRoT-7 system are reported in Table |4] Considering our value 
of the mass of CoRoT-7b, one sigma higher than the published 
value, and its refined measurement of radius, slightly smaller 
= 1.58 + 0.10 REarth, derived by Bruntt et al. (2010) using im- 
proved stellar parameters, Valencia et al. (2010) predicted that is 
a planet compatible with an Earth-Uke composition (33% iron, 
66% silicate). 

3.4. L Hor 

iHor, or HD 17051, is a young GOV star with V = 5.40. The 
320. 1-d period planet tHorb with a small eccentricity of 0.161 
was reported by Kiirster et al. (2000). They noted an excess RV 
scatter of 27 ms"' due to stellar activity. Rocha-Pinto & Maciel 
(1998) measured a Call index of logR'^^--4.65, and Jenkins et 
al. (2006) measured log7?^^=-4.59. Asteroseismologic observa- 
tions were made with the high-precision spectrograph HARPS 
of I Hor (Vauclair et al. 2008). They studied the acoustic oscil- 
lations of the star and demonstrated that i Hor was formed in 
the same primordial cloud as the Hyades. The data helped re- 
fine the mass of the star M» = 1.25 + 0.01 M©, which lead us to 
re-evaluate the i Horb planetary mass to be 2.6 M jup. 

We studied these data to characterize the active jitter and 
search for a possible hidden Doppler motion. Observations were 
carried out during eight consecutive nights between 19 and 26 
November 2006 with 1856 measurements each of time exposure 
of about 100s. The p-modes are visible at high frequencies and 
a low frequency RV signal with amplitude of about 20 ms"' ap- 
pears because of the activity. We kept 1659 measurements with 
SNR(550nm) > 120. We averaged the data in groups of 20 mea- 
surements in order to average the p-mode signature. The mean 
RV photon-noise uncertainty in the averaged points is then about 
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Fig. 19. RVhigh and RViok of ( Hor derived from HARPS spectra 
as a function of time. The slope due to the known planet has not 
been removed. The error bars are also plotted. 



26 cms ' , but the true precision is limited by the instrumental ac- 
curacy X! 80 cms"'. 

Before studying the RV variations due to stellar activity, we 
subtracted the long -period planet Doppler motion. We computed 
the L Hor b Keplerian solution with the parameters of Kiirster et 
al. (2000) because our data has a too short duration and cannot 
be used to better constrain the planetary orbit. We subtracted 
a slope of -650 ms 'yr ' from the RV measurements. The 
uncertainty in the Kiirster et al. (2000) parameters, particularly 
in the period {SP = 2. Id) and the epoch of maximum RV (dTo = 
3.0 d), induced an uncertainty of +50 ms^'yr"' in the value of 
the slope at the epoch of our observations, approximately eight 
years after the Kurster ones. 

To check that the main RV variations is due to activity, we 
computed and compared RVhigh and RViow in Fig- [IS] The am- 
plitude of RWiow is roughly 1.5 times the amplitude of RVhigh, 
as expected in the case of stellar activity (cf. Sect. 12.31 1. Thus, in 
Fig.|20]an anti-correlation is observed between the RV and the 
'Vspan in agreement with an active signature. 

The study of the dark spot simulations in Sect. |2] shows that 
the active RV jitter is well-fitted when the rotational period of 
the star is known. The t Hor rotational period was estimated by 
Saar & Osten (1997) and by Saar et al. (1997) from Call emis- 
sion. They found Proi - 7.9d and 8.6d, respectively. From the 
Kurucz model atmosphere, Kurster et al. (2000) derived a vsin/ 
- 5.5 kms"', which agrees with the 5.7 kms"' found by Saar 
& Osten (1997). Also in agreement, we derived from the mean 
FWHM of the HARPS CCF a value of v sin / = 5.7 kms '. Using 
the Saar & Osten (1997) formula to estimate the sin / from the 
V sin i, Prot, and the stellar radius, we derived the stellar inclina- 
tion /=55+10°. Following Saar & Donahue (1997), we estimate 
that the /j, nonuniform portion of the spot distribution responsi- 
ble for the observed variability is fs 

We choose arbitrarily a rotational period equal to 8.2 days. 
The RV are fitted by three sinusoids with periods fixed at the 
rotational period and its two first harmonics (cf. Fig. ISTT l. The 
residuals are equal to cr = 1.03ms"' reaching almost the instru- 
mental accuracy. A peak at the third harmonic is barely detected 
above the noise in the Lomb-Scargle periodogram of the residu- 
als, as expected when data still exhibit signature of activity (cf. 
Sect. 12.5b . Fitting simultaneously the fourth sinusoids, the semi- 
amplitude of the fourth one is only 1.5 ms"' and the residuals 
are slightly smaller at a value of cr =0.87ms"'. The amplitude of 
this last sinusoid represents the limit of our detection (< 2 cr) as 
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rived from HARPS spectra. The line is the least squares fit. The 
ranges have the same extents along the x- and y-axes. The er- 
ror bars are also plotted. One may compare this shape with that 
of the simulation of two spots separated by 120° in longitude in 
Vm.\n[ibottoin). 
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Fig. 21. Top: RV of i Hor derived from HARPS spectra as a 
function of time. The three-sinusoid fit, with periods fixed at the 
rotational period of the star and its two first harmonics, is plotted 
as a solid line. The rotational period is chosen to be equal to 8.2 
days. Bottom: Residuals from the fit as a function of time. The 
dispersion is equal to 1.03 ms"'. 



before in the study of the dark spot simulations, illustrating that 
three sinusoids is sufficient to remove most of the active jitter. In 
Fig. |22] it is shown that minimum values of are obtained for 
Prot between 7.9d and 8.4d. 

In the iHor data, we do not detect a short-period companion. 
Nevertheless, we wished to check whether we have subtracted 
the RV shift due to a companion when subtracting the effects of 
activity. We ran simulations and added RV due to fake planets 
to the tHor data. We consider only the case of circular orbits, 
which is a good hypothesis for planets of period shorter than 
6 days (Eggenberger & Udry 2010). We fitted the active jitter 
with three sinusoids of period fixed to the rotational period and 
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Fig. 22. Reduced as a function of the value of the rotational 
period used for the fit of the ( Hor data. The rotational period is 
used to fix the period of the fitted sinusoids. Minimum values of 
the;^^ are obtained for Prot between 7.9d and 8.4d. 
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Fig. 23. Minimal semi-amplitude of a planet that would be de- 
tected in the iHor data as a function of the planetary period. The 
error bars represent the difference between the fitted and the sim- 
ulated parameters, or, if it is greater the error in the fitted param- 
eters. 



its two first harmonics. Afterwards, we considered the Lomb- 
Scargle periodogram of the residuals. If a peak at the planetary 
period was detected, we fitted simultaneously a Keplerian with 
null-eccentricity to obtain the planetary parameters. We consider 
that a peak is significant if its false-alarm probability is smaller 
than 10"^. In Fig.|23] we plotted the minimal semi-amplitude de- 
tected in the iHor data as a function of the planetary period. The 
error bars represent the difference between the planetary period 
and semi-amplitude that were fitted and the simulated ones or, if 
it is greater, the error in the fitted parameters. Periods longer than 
7 days were not considered because observations during at least 
one planetary period were needed to detect a planet. Planetary 
periods between 7 and 2.5 days are not easily characterized be- 
cause, without other high-precision measurements, we cannot 
evaluate the amplitude of the RV jitter due to activity and then, 
estimate above which amplitude the RV variation might be due 
to a companion. When shorter than 2.5 days, the periods are de- 
tected in the periodogram of the residuals, after subtraction of 
the RV jitter. 

To conclude, the (Hor RV modulations can be clearly ex- 
plained by dark spots on the photosphere of a star with a rota- 
tional period in the range [7.9-8.4] days. There is no doubt that 
the main RV variations are due to the stellar activity by the anti- 



I. Boisse et al.: Disentangling between stellar activity and planetary signals 



11 



correlation RV-V.5y,n„ (Fig. l20b and the difference in amplitude 
between the RV/ov, and RV /„;,/, (Fig. \T9[ . We excluded the pres- 
ence of planets with minimum masses between 6 and 10 MEarth 
and respective periods between 0.7 and 2.4 days (Fig. l23T l. The 
causes limiting our diagnostic are primarily the uncertainty in 
the rotational period of the star, the short observational period, 
and, to a lesser extent, the uncertainty in the slope due to the 
iHorb planet. 

4. Conclusion and Perspectives 

We have performed simulations of stellar dark spots to check 
the validity of the standard diagnostics of activity and to derive a 
method to subtract RV variations due to activity. We have stud- 
ied the effects of stellar spots on photometry, measured radial 
velocities, and other spectral parameters, including the asymme- 
try of the mean line of the spectrum that we estimate with a new 
parameter, span-^^ high - RV/ow, which is less sensitive to noise 
than the usual bisector span BIS. We have shown that the anti- 
correlation between RV and the Vjpo,, needs a minimal sampling 
of 0.1 of a stellar rotational period to keep its diagnostic power 
for several spots, and that an anti-correlation might be observed 
between the photometry and the FWHM of the spectral lines. We 
have pointed out that bright spots also induce an anti-correlation 
between RV and the V span- We have shown that stellar spots in- 
duce RV variations at the rotational period Prot of the star and its 
two first harmonics Protl'^ and Protl^- These RV variations can 
be corrected by three sinusoids, the remaining activity signal at 
the following harmonic representing less than 10 % of the initial 
RV amplitude. Fitting stellar activity with three sinusoids may 
be used to disentangle it from Doppler motion provided that 1) 
the period of the planet is not close to that of the stellar rota- 
tion or one of its first harmonics, and 2) the rotational period of 
the star is accurately known or the data cover at least two ro- 
tational periods of the star. Under these conditions, this method 
could even reveal low-mass planets with semi-amplitudes down 
to about 30% of the semi-amplitude of the activity-induced vari- 
ation. 

These results have been validated on four known active 
planet-host stars. In the case of CoRoT-7, we have also detected 
the two planets reported by Queloz et al. (2009) but our uncer- 
tainties in the fitted parameters are several times larger to prop- 
erly take into account the stellar active jitter Moreover, with our 
simultaneous modeling of the activity and planetary parameters, 
slightly higher masses are found: 5.7 + 2.5 MEanh for CoRoT-7b 
and 13.1 +4.1 MEarih forCoRoT-7c (assuming that both planets 
are on coplanar orbits). In the case of t Hor, we have been able to 
exclude low-mass planets with periods between 0.7 and 2.4 days 
with semi-amplitude greater than 4 ms"', which correspond re- 
spectively to a minimum mass of between 6 and 10 MEanh- 

Our simulations did not account for the differential rotation. 
This may generate frequencies other than the harmonics of the 
rotational period. These additional periodicities, however, were 
not detected in our study of known active stars. On the other 
hand, for an RV survey, when a periodicity close to the rota- 
tional period of the star or its harmonics is found, a photometric 
follow-up has to be performed. As illustrated in the discovery 
and confirmation of the planetary companion of HD 192263, a 
planet may be found with a period roughly close to that of the 
photometric variations of the star (Santos et al. 2000; Henry et 
al. 2002; Santos et al. 2003). Simultaneous photometry with the 
RV is then needed to check the phase of the parameters. 

This method for distinguishing stellar activity and planetary 
signals is particularly suitable for RV follow-up of transit sur- 



veys from the ground or space such as CoRoT or Kepler In this 
case, photometry of active stars is indeed well sampled and char- 
acterized (Paulson et al. 2004). The method could also be used 
for RV survey around (low) active stars when the companion 
candidate has a periodicity that differs from the rotational period 
of the star and its harmonics, to refine the planetary parameters. 
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Table 1. Sinusoidal orbital solutions for fitting RV jitter activity for different inclinations of the star and latitudes of one stellar spot 
(Fig.[B. 



Parameters 


i=90° and lat= 


=+30° 


Pp 
K [ms-'l 
o"(o-c) [ms"'] 


Pro, 

24.83 ± 0.02 


Pro,/2 

32.65 ± 0.02 
before the fit : 33.77 


Pro,/3 

22.78 ± 0.02 
after the fit : 1.26;^^ = 


1.22 


Pro,/4 

8.66 ± 0.02 


Pro,/5 

1.14 + 0.03 


1=40° and lat= 


=+30° 


Pp 
K [ms-'l 
(Tio-o [ms-'] 


30.97 ± 0.02 


Pro,/2 

26.68 ± 0.02 
before the fit : 29.58 


Pro,n 

8.57 ± 0.02 
after the tit : 1.05;^^ = 


1.01 


Pro,/4 

1.79 ±0.04 


Pro,/5 

0.85 ± 0.03 


i=40° and lat= 


=+50° 


Pp 
K [ms-'l 
(T(o-o [ms-'l 


Pro, 

40.59 ± 0.02 


Pro,/2 

20.80 ± 0.02 
before the tit : 32.41 


Pro,n 

4.29 ± 0.02 
after the tit : 1.04;^^ = 


1.00 


Pro,/4 

0.97 ± 0.02 


Pro,/5 

0.18+0.02 



Table 3. CoRoT-7 fit with three sinusoids for the active jitter and two Keplerians for the planets for the last 37 days of HARPS RV 
measurements. 



Parameters 



Pp [days] 


Pro,=23 (fixed) P,„, 12= 1 1 .5 (fixed) 


P„, 13=1 .66661 (fixed) 


3.695 + 0.02 


0.8536 (fixed) 


K[ms-'] 


14.3 + 0.8 3.7 + 0.8 


1.2 + 0.2 


6.1 +0.6 


4.5 + 0.7 


e 






(fixed) 


(fixed) 


To [JD] 






54899.2 + 0.7 


54899.761 (fixed) 


o-(o-c) [ms-'] 


before the fit 10. 


after the fit 3.1; reduced;^^^ = 2. 







